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Abstract 

We present a brief introduction to the physics of parton saturation/Color 
^ I Glass Condensate (CGC). 
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^ . 1. Introduction 

One of the most interesting features of quantum chromodynamics (QCD) 
is the property of the asymptotic freedom: the strong coupling constant 

»vq ' is small at large momenta/short distances, and it is large at small mo- 

O ■ menta/large distances |l|, |2|. Finding the scale that determines the value of 

the characteristic running QCD coupling is one of the the central questions 
for high energy scattering physics, important for the theoretical description 

O ! of both the hadronic and nuclear scattering processes. 

A naive answer to this question would be to say that in high energy 
scattering the large center-of-mass energy s determines the scale of strong 

^ . coupling making it small: as{s) <C 1. While such statement would be true for 

several s-channel processes, the dominant contribution to total cross sections 
in high energy hadronic and nuclear scattering comes from the t-channel ex- 
changes, for which the scale of the QCD coupling constant is not given by s, 
but by the typical transverse momentum in the problem. For proton-proton 
scattering one may estimate the typical transverse momentum to be of the 
order of the inverse transverse size of the protons, which is roughly the QCD 
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confinement scale Aqcd- People performing such estimate would pessimisti- 
cally conclude that the QCD coupling in high energy hadronic scattering runs 
as asiAg^^) and is therefore not small, as{AQQ^) ~ 1. With the coupling 
constant of order-one, we would have little chance of describing the total 
hadronic and nuclear cross sections from first principles, at least certainly 
not with the help of QCD perturbation theory. 

Traditional perturbative QCD (pQCD) approaches are well-aware of the 
above problem, and try to avoid it by separating hard sub-events where 
the coupling is small from the full event with large QCD coupling. For in- 
stance, in hadronic scattering pQCD may be used to calculate jet production 
cross section, where the hard partonic scattering is factorized from the non- 
perturbative distribution and fragmentation functions. The high transverse 
momentum pt of the produced hard parton insures applicability of pQCD: 
C(s{Pt) ^ 1- Similarly, in deep inelastic scattering (DIS) pQCD can describe 
structure functions at high photon virtuality Q^, since there as{Q^) <C 1, but 
pQCD is expected to fail at low-Q^. Since jet production events in hadronic 
collisions constitute a small percentage of the total cross section, pQCD ap- 
proach describes only rare events, and is not applicable to the description of 
the bulk particle production and dynamics. 

Saturation physics provides a new way of tackling the problem of total 
hadronic and nuclear cross sections. It is based on the theoretical observation 
that small-a; hadronic and nuclear wave functions, and, therefore, the scat- 
tering cross sections as well, are described by an internal momentum scale 
known as the saturation scale and denoted by Qs [3j] • This intrinsic momen- 
tum scale grows with the center-of-mass energy s in the problem, and with 
the increasing atomic number of a nucleus A (in the case of a nuclear wave 
function) approximately as 

Ql ~ A'/'s' (1) 

where the best current theoretical estimates of A give A = 0.2 -^ 0.3 j^. 
Therefore, for hadronic collisions at high energy and/or for collisions of large 
ultrarelativistic nuclei, saturation scale becomes large, Q^ ^ Aqq^. Since 
for total cross sections Qs is usually the only momentum scale in the problem, 
we expect it to give the scale of the running QCD coupling constant, making 
it small 

«.(Q') « 1 (2) 



and allowing for first-principles calculations of total hadronic and nuclear 
cross sections, along with extending our ability to calculate particle produc- 
tion and to describe diffraction in a small-coupling framework. 

Below we present a short review of the saturation physics, which is also 
known as the Color Glass Condensate (CGC) physics. For more extensive 
descriptions of the subject we refer the readers to the review articles |5|, la, LZJ- 

2. Brief Review of Saturation Physics 

Traditional approach to saturation physics consists of two stages, corre- 
sponding to two different levels of approximation. The first approximation 
corresponds to the classical gluon field description of nuclear wave functions 
and scattering cross sections. It resums all multiple rescatterings in the nu- 
cleus, but lacks rapidity-dependence. The latter is included through quantum 
corrections, which are resummed by the non-linear evolution equations. This 
constitutes the second level of approximation. We will present both stages 
below. 

2.1. Classical Gluon Fields 

Imagine a single large nucleus, which was boosted to some ultrarelativis- 
tic velocity, as shown in Fig. [TJ We are interested in the dynamics of small-x 
gluons in the wave function of this relativistic nucleus. The small-x glu- 
ons interact with the whole nucleus coherently in the longitudinal direction: 
therefore, only the transverse plane distribution of nucleons is important for 




Figure 1: Large nucleus before and after an ultrarelativistic boost, 
the small-x wave function. As one can see from Fig. [1], after the boost the 



nucleons, as "seen" by the small-x gluons, appear to overlap with each other 
in the transverse plane, leading to high parton density. Large occupation 
number of color charges (partons) leads to classical gluon field dominating 
the small- a; wave function of the nucleus. This is the essence of the McLerran- 
Venugopalan (MV) model [8|. According to the MV model, the dominant 
gluon field is given by the solution of the classical Yang-Mills equations 

v^F^" = r (3) 

where the classical color current J'^ is generated by the valence quarks in the 
nucleons of the nucleus from Fig. [1] 

The equations ([3]) were solved for a single nucleus exactly |9|, llO| result- 
ing in the unintegrated gluon distribution 0(x, k"^) (multiphed by the phase 
space factor of the gluon's transverse momentum kx) shown in Fig. |2] as 
a function of kx- (Note that in the MV model 0(x, k"^) is independent of 
Bjorken-x.) Fig. [2] demonstrates the emergence of the saturation scale Qs'- 
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Figure 2: Unintegrated gluon distribution (j){x, k^) of a large nucleus due to classical gluon 
fields (solid line). Dashed curve denotes the lowest-order perturbative result. 

as one can see from Fig. [21 the majority of gluons in this classical distribution 
have transverse momentum kx ^ Qs- Since in this classical approximation 
Ql ~ A^/^, for large enough nucleus all of its small-x gluons would have 
large transverse momenta kx ^ Qs ^ ^qcd, justifying applicability of per- 
turbative approach to the problem. Note that the gluon distribution slows 
down its growth with decreasing kx for kr < Qs (from power-law of kx to a 
logarithm): the distribution saturates, justifying the name of the saturation 
scale. 



2.2. Nonlinear Evolution 

While the classical gluon fields of the MV model exhibit many correct 
qualitative features of saturation physics, and give predictions about A- 
dependence of observables which may be compared to the data, they do 
not lead to any rapidity/ Bjorken-x dependence of the corresponding observ- 
ables, which is essential in the data on nuclear and hadronic collisions. To 
include rapidity dependence one has to calculate quantum corrections to the 
classical fields described above. 

new parton is emitted as energy increases 
it could be emitted off anyone of the N partons 



N partons any two partons can recombine into one 



Figure 3: Nonlinear small-x evolution of a hadronic or nuclear wave functions. All partons 
(quarks and gluons) are denoted by straight solid lines for simplicity. 

The inclusion of quantum corrections is accomplished by the small-x evo- 
lution equations. The first small- a; evolution equation was constructed be- 
fore the birth of the saturation physics. This is the Balitsky-Fadin-Kuraev- 



Lipatov (BFKL) evolution equation ll|, |l2|. This is a linear evolution equa- 



tion, which is illustrated by the first term on the right hand side of Fig. [31 
Consider a wave function of a high-energy nucleus or hadrons: it contains 
many partons, as shown on the left of Fig. O As we make one step of evo- 
lution by boosting the nucleus/hadron to higher energy, either one of the 
partons can split into two partons, leading to an increase in the number of 
partons proportional to the number of partons A^ at the previous step, 

dNix,^^) _ ^ ^^ 



91n(l/x) 



Oi,KBFKL®N{xXT). (4) 



with Kbfkl an integral kernel. Clearly the BFKL equation (jl]) introduces 
Bjorken-x/rapidity dependence in the observables it describes. 

The main problem with the BFKL evolution is that it leads to the power- 
law growth of the total cross sections with energy, atot ~ s"^~^, with the 



BFKL poineron intercept ap — 1 = {AagNc ln2)/7r > 0. Such power-law 
cross section increase violates the Froissart bound, which states that the total 
hadronic cross section can not grow faster than In s at very high energies. 
Moreover, power-law growth of cross sections with with energy violates the 
black disk limit known from quantum mechanics: high-energy total scattering 
cross section of a particle on a sphere of radius R is bounded by 

atot <2nR\ (5) 

(Note the factor of 2 which is due to quantum mechanics, this is not simply 
a hard sphere from classical mechanics!) 

We see that something has to modify Eq. (^ at high energy. The modi- 
fication is illustrated on the far right of Fig. [31 at very high energies partons 
may start to recombine with each other on top of the splitting. The recom- 
bination of two partons into one is proportional to the number of pairs of 
partons, which, in turn, scales as A^^. We end up with the following non- 
linear evolution equation: 



dN{x,k 



91n(l/x) 



as Kbfkl ® N{x, kl) - as [N{x, 4)]^. (6) 



This is the Balitsky-Kovchegov (BK) evolution equation |13l . Il4j |. which is 
valid for QCD in the limit of large number of colors A''^. An equation of 
this type was originally suggested by Gribov, Levin and Ryskin in |3| and 



by Mueller and Qiu in [15[, though at the time it was assumed that the 



quadratic term is only the first non-linear correction with higher order terms 



possibly appearing as well: in |l3l . |14| the exact form of the equation was 



found, and it was shown that in the large- A"c limit Eq. (jH]) does not have any 
higher-order terms in A^. Generalization of Eq. (^ beyond the large- A^c limit 
is accomplished by the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov- 
Kovner ( JIMWLK) |16l . llTj evolution equation, which is a functional differ- 



ential equation. 

The physical impact of the quadratic term on the right of Eq. (E]) is clear: 
it slows down the small-x evolution, leading to parton saturation and to total 
cross sections adhering to the black disk limit of Eq. (|5]). The effect of gluon 
mergers becomes important when the quadratic term in Eq. ([6]) becomes 
comparable to the linear term on the right- hand- side. This gives rise to the 
saturation scale Qs, which now grows with energy (on top of its increase with 
A), as was advertised in Eq. ([T]) above. 
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Figure 4: Map of high energy QCD in the (Q^, Y = ln 1/x) plane. 



We summarize our knowledge of high energy QCD in Fig. |H in which 
different regimes are plotted in the (Q^, 1^ = In 1/x) plane, by analogy with 
DIS. For hadronic and nuclear collisions one can think of typical transverse 
momentum p^ of the produced particles instead of Q^. Also rapidity Y and 
Bjorken-x variable are interchangeable. On the left of Fig. IHwe see the region 



with Q < Aq(jjj in which the coupling is large, 



a.. 



1, and small-coupling 



approaches do not work. In the pessimistic view of high energy scattering 
described in the Introduction, this is exactly where the total hadronic and 
nuclear cross sections would be. In the perturbative region, Q^ ^ ^qcd^ 
we see the standard DGLAP evolution and the linear BFKL evolution. The 
BFKL equation evolves gluon distribution toward small-a;, where parton den- 
sity becomes large and parton saturation sets in. Transition to saturation 
is described by the non-linear BK and JIMWLK evolution equations. Most 
importantly this transition happens at Ql ^ ^qcd where the small- coupling 
approach is valid. 

2.3. Some CGC Phenomenology 

One of the important predictions of saturation/CGC physics was the so- 
called geometric scaling of the total DIS cross section. It was argued that 



with Qs{x) being the only scale in the problem, DIS structure functions and 
cross sections should depend only on one variable - Q'^/Ql{x), instead of 
being functions of two variables x and Q^. This prediction was supported by 
detailed calculations based on BK evolution [l8|, [l9| and was confirmed by 



an analysis of the HERA DIS data (20 
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Figure 5: Nuclear modification factor as a function of kx/Qs as predicted by satura- 
tion/CGC physics. Different curves correspond to different rapidities, with lower curves 
corresponding to higher rapidity. 

Another prediction of the non-linear evolution (jH]) concerns particle pro- 
duction in proton-nucleus (d+Au) collisions. It follows from Eq. ([6]) combined 
with the formulas for particle production in CGC that the nuclear modifi- 
cation factor RP^ should decrease as one goes toward more forward rapidity 
2lL I22L |23| . This prediction is illustrated in Fig. |5] and was confirmed by 



RHIC experiments. 



3. Recent Progress in Saturation Physics 

In recent years saturation physics has become more precision-oriented, 

with steps taken to improve its quantitative predictions. Running coupling 

corrections to the BFKL/BK/JIMWLK evolution have been calculated in 

241 |25| . It led to an interesting result, in which the fixed coupling as in 

Eq. (|6]) was replaced by a "triumvirate" of the running couplings at different 



scales 24 l25 : 



Such behavior has never been seen in field theories outside small- a; physics. 
Among other recent developments, next-to-leading-order BK equation has 



been found in [26[, resulting in a rather complicated but useful expression. 

There is hope that we have finally managed to significantly advance both 
the qualitative and quantitative understanding of QCD at high energies. 
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